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Abstract: We consider the creation of the cosmological baryon asymmetry in the Two 
Higgs Doublet Model. We imagine a situation where the masses of the five Higgs particles 
and the two Higgs vevs are constrained by collider experiments, and demonstrate how the 
requirement of successful baryogenesis can be used to further constrain the remaining 4- 
dimensional parameter space of the model. We numerically compute the asymmetry within 
the scenario of Cold Electroweak Baryogenesis, which is particularly straightforward to 
simulate reliably. 
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1 Introduction 

With the discovery of what is expected to be the Higgs particle [2], the Standard Model of 
particle physics now provides a coherent and consistent theory of fundamental physics up to 
and including the electroweak scale. Although many phenomena in the realm of cosmology, 
such as inflation, dark matter and dark energy, are not addressed in this framework, all 
experimentally observed processes are very well described 1 . 

One central issue at the boundary between cosmology and particle physics is the origin 
of the baryon asymmetry observed in the Universe. A substantial effort has been made to 
link this phenomenon to electroweak scale physics [3-5] , since it is the lowest energy where 
baryon number violation may occur. In combination with C-, P- and CP-violation and 
an out-of-equilibrium electroweak symmetry breaking transition, baryon asymmetry with 
the observed magnitude can indeed be produced. However, because Standard Model CP- 
violation is minute [6] (see however [7]) and the electroweak transition is a crossover in the 
Standard Model for the physical Higgs mass of 125-126 GeV [8], electroweak baryogenesis 
requires additional fields and interactions to exist. 

1 We will consider the Standard Model to include right-handed neutrinos and a non-zero mass term for 
neutrinos. Hence the phenomenon of neutrino oscillations is considered part of the Standard Model. 
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Presumably the simplest way to achieve successful electroweak baryogenesis is to ex- 
tend the scalar sector of the Standard Model by an additional field. This could be a SU (2) 
singlet but probably the most popular extension is the Two-Higgs Doublet Model (2HDM) 
with an additional SU{2) doublet. A number of different "types" exist, depending on how 
the two Higgs doublets couple to the fermions (see [9] for a recent review) . The most general 
Higgs potential contains 14 (real) parameters, including up to two CP violating phases. 
In addition, because of the richness of the vacuum structure CP may be spontaneously 
broken. 

Restricting to a sub-class of models with only 10 real parameters, we ask the question 
whether the observed baryon asymmetry can be used to constrain the parameter space, 
complementing direct collider experiments. We expect that masses (4 different, of which 
we know the lightest) and vevs (2, of which we know one) are the easiest to measure, and so 
we will imagine that in future these are constrained, leaving a 4-dimensional less accessible 
subspace. Our aim here is to show how one may in principle sweep through this subspace 
and potentially use the observed baryon asymmetry to pin down the allowed parameter 
region. 

The 2HDM doublet can accommodate a strong first order phase transition, but we 
will consider a different scenario, where electroweak symmetry breaking is a cold spinodal 
transition [10-13]. This is a viable alternative to the standard "Hot" scenario [4], but for 
our purpose here, its main virtue is that it is straightforward to simulate numerically from 
first principles. Cold electroweak baryogenesis may be realized as a result of coupling to 
another scalar field, which may [14] or may not [15, 16] be the inflaton. We demonstrated 
in [1] through direct numerical simulations that a baryon asymmetry is indeed produced, 
as a result of the interplay between the explicit CP/C violation in the Higgs potential 
and the C- and P-violating (but CP-conserving) gauge- fermion interactions. It turns out 
that when simulating a bosonized version of the theory, it is necessary to include the P- 
breaking of fermions, and we did this through an effective higher order bosonic interaction 
term, parameterized by a coefficient #c/P- Although this coefficient can in principle be 
computed analytically, this is a very non-trivial task and we chose to keep it as a free 
parameter. We found that in order for the observed asymmetry to be reproduced, we need 
6c/p — (2 to 3) x 10 -4 or larger. 

The paper is structured as follows: In section 2 is a brief introduction of the bosonized 
2HDM. In section 3 we parametrize the 4-dimensional parameter space in terms of a field 
transformation, two angles and one mass scale. The numerical results are presented in 
section 4, where we for a given range of parameters, and using lattice simulations in real- 
time, directly compute the baryon asymmetry in the bosonized electroweak sector. We 
conclude in section 5. In Appendix A, we further discuss the parametrization of the neutral 
Higgs masses. 
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2 The 2HDM 



The 2HDM is denned through the continuum action 
S = - I d 3 xdt 



-l^TrF^F^ + (D^Dffa + {D^D^fa + V(4> u 2 ) + C c/P 



(2.1) 



where we use the metric rj = diag( — \- ++), 0i 2 are SU(2) doublets with hypercharge +1 
and F^ u is the field strength tensor of the gauge field. We will ignore the SU(3) and U(l) 
gauge fields. The covariant derivative is D^fa = (9 M + iA^fa and the potential is in all 
generality 2 
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(2.2) 

The parameters Ai >2 ,3,4 and //f 1 22 are real and in general ^5^,7 and // 2 2 are complex. In this 
paper, we only study the 2HDM with a softly broken Z 2 symmetry, in which Xq = A7 = 
[9]. There is then only one independent CP violating phase. In the Standard Model as well 
as the 2HDM there is also CP-violation through the complex phase in the CKM mixing 
matrix. For the purpose of the present work, we will assume that the effective CP-breaking 
terms arising from this are negligible, although at very low temperatures, this may not be 
correct [7]. 

We will take (2.1) to represent a bosonized version of the full theory, where fermions 
have been integrated out, and their effect is captured in a C- and P-breaking term given 

by [17] 

l w 

The Yukawa couplings and the mixing matrix is encoded in the real parameter 5q/p, and 
it can in principle be computed from the model. The standard prescription in bosonized 
theories, which we will also adopt here, is then to infer the value of the baryon number B 
through the anomaly equation 

B(t) - B(0) = 3[N cs (t) - N cs (0)], (2.4) 

where N cs is the Chern-Simons number of the SU(2) gauge field. 

The reason for including the term (2.3) is that, as demonstrated in [1], to generated 
a non-zero average Chern-Simons number, we need P-symmetry to be broken as well as 
CP-symmetry. It is easy to see that (2.3) conserves CP. 



2 We here correct an error in [1] in the normalization of the coefficients. The results obtained there were 
based on the conventions presented here. 
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1 2 

It turns out that in a finite temperature environment, the Higgs winding numbers N w ' 
for the two Higgs fields, respectively, are much cleaner observables. At late times, the three 

1 2 

agree, N w ' = N cs , and so we will identify the winding numbers at the end of the simulation 
to be the late time value for Chern-Simons number and hence the baryon asymmetry. 

3 Choices of the parameters 

3.1 The full parameter space of the 2HDM 

We will re-parametrize the 10-dimensional parameter space in the following way: 

• Vacuum parameters (3): v, (3 and 6. 

Without loss of generality, we can parametrize the Higgs fields in terms of 2 complex 
and 4 real degrees or freedom as 




<h = 

vi +m +«Xi)/v / 2 / ' \ (v 2 + r/2 + /\/2 



(3.1) 



and 



Xi = cos/3G° -sin/37? 3 , X 2 = cos/3r/ 3 + sin^G , (3.2) 
= C os (3G + - sin (3H + , 0+ = cos f3H + + sin f3G + . (3.3) 

The vacuum is given by = (j}f 2 = ? 7i,2,3 = 0, in terms of v\e td and V2- We 

introduce v and /3 through 

i>i = tjcos/3, V2 = vsinf3, v<ilv\ = tan/3. (3.4) 

Minimizing the Higgs potential gives three equations 

9 d d 

-q^-V\vi,V2,9 = 0) ^~MW1,U2,0 = 0' ^^1^1,112,6=0, (3.5) 

with which we can replace three couplings/mass parameters by /3, 6 and v. 
Higgs masses (4) : mi 2,3 and m±. 

There are five physical Higgs bosons: two form one charged field and the rest 
are mass eigenstates formed as linear combinations of the neutral fields 771,2,3- We 
introduce the mass eigenvalues for these, m± and 77112,3, respectively, and these 
replace four other parameters (see also Appendix A). 

Neutral Higgs mixing angles (2) : ax, oli- 

As discussed in Appendix A, the mass matrix of the neutral Higgs modes is in general 
not diagonal in the fields 7/1,2,31 but it can be diagonalized through three mixing 
angles 0^2,3- Only two of these are independent, and we take 03 to be fixed through 
Eq. (A. 10), which has 0, 1 or 2 solutions for a given set of («i,a2). 

A mass parameter (1) : fi 2 = Re( / uf 2 e_i61 )- 
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At the end of the day, the parameter set in the Higgs potential denoted by {A} is a 
function of the above 10 parameters. In the following discussion, for simplicity of notation 
we will use {A}[. . .] with the ellipsis being some of the above parameters relevant for the 
discussions only. 

3.2 The subspace spanned by (ai, a.2, 6, fi) 

As explained in the Introduction, we will assume that the 4 distinct masses and the two 
vevs have been determined (or at least constrained) by experiment, so that we can assign 
values to them: 

• Vacuum parameters: 

The vev v is known but not j3 and we choose 

u = 246 GeV, tan/3 = 2. (3.6) 

• Higgs masses: The lowest neutral Higgs mass mi is fixed by experiment [2]. Based 
on our choice of /3 and experimental constraints [9], we choose 

mi = 125 GeV, m 2 = 300 GeV, m 3 = 350 GeV, m± = 400 GeV. (3.7) 

This leaves a 4-dimensional parameter space, spanned by a±, a 2 , 6, and fi. 
3.2.1 Symmetries 

Symmetries in the Higgs potential V help us to further simplify our calculations. Since V 
is real, it follows that 

{A}^{A}*, Si-^SJ, (3.8) 
is a symmetry. This imposes the relation between the sets of parameters {A}, 

{X}[ai, a 2 ,a 3 ] = ({A}[«i, -a 2 ,vr - a 3 ])* , (3.9) 

which is equivalent to the charge conjugation of the bosonic fields according to (3.8). There- 
fore, the generated baryon asymmetry flips sign when complex conjugating the parameter 
set {A}, i.e., 

N cs [ai,a 2 ,a 3 ] = -N cs [a>i,-a 2 ,Tr - a 3 ], 

Nl' 2 [a 1 ,a 2 ,a 3 ] = -N^ 2 [ati, -a 2 , 7r - a 3 ], (3.10) 

and so there is a redundancy between the upper and lower half-plane in a\ — a 2 space. 
Finally, the symmetry 

fa^e-'Ui, X 5 ^e~ 2ie X 5 , mL^^Vi2- (3-H) 

will also be very useful. Using this transformation, one can easily see that 

A 5 [0]=e^A 5 [O], M? 2 [^=e i V? 2 [0]- (3.12) 
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Therefore, one can first find the parameter set {A}[0] and then obtain {A}[#] by the above 
transformation. The physical Higgs masses are unchanged under such a transformation. 
Since fi is invariant under (3.12), we need only consider varying the potential in the 3- 
dimensional parameter space spanned by (a±,a2,fJ-), and we get the ^-direction for free. 

The potential at different values of 9 are equivalent, but with different field basis. 
However, the symmetry in (3.12) is explicitly broken as soon as the scalar sector is coupled 
to fermions, or in our case the C-/P-violating term in (2.3) is included. Then different 8 
are physically distinct, as under the transformation (3.11), 



3.2.2 Basic constraints, maxima and saddle points 

For a given value of /i, we now survey the whole a\ — oti plane. For each such pair, we 
accept/reject based on overall stability (potential is bounded from below), unitarity (tree- 
level Higgs-Higgs scattering amplitudes are smaller than unity), and whether the minimum 
found is a global minimum. Conditions for stability and unitarity are well-known and the 
interested reader is referred to [9] and references therein. For the requirement of the global 
minimum, we find all the other minima of the potential and establish that the chosen one 
is in fact the one with lowest potential energy. We also reject if there are no solutions for 
Ck3, and finally we reject if the potential has a minimum at (</»i, 02 ) = (0, 0) (see below). 

The origin (0i, 02) = (0, 0) is always a stationary point of V. For each of the surviving 
pairs (ai,a2), we compute the eigenvalues of the mass matrix at the origin (not to be 
confused with M 2 of (A. 7), the neutral Higgs sector mass matrix in the minimum), 



If both eigenvalues are negative, both Higgs fields will experience a spinodal transition, 
and we name this parameter point a maximum. If only one eigenvalue is negative (and the 
other positive), only one field goes spinodal, and we name the parameter point a saddle 



point. In principle, such a minimum could lead to tunneling and bubble nucleation on 
the way to symmetry breaking, but this returns us to standard electroweak baryogenesis, 
which we do not consider here (but see also [16]). 

Fig. 1 shows the a\ — ct2 plane for various values of fi, where we have sampled points 
with a spacing of tt/40. We have indicated maxima by fat black dots, and saddle points 
by smaller red dots. The rest of the parameter space has been discarded for one of the 
reasons explained above. Where a red and a black dot are superposed, this corresponds to 
the two values of Q3, and that these give a maximum and a saddle point, respectively. The 
lines 02 = 0, ±7r/2 have zero CP-violation, and can therefore not provide baryogenesis. 

At zero fi, no choice of a± t 2 survives the constraints. For small, but non-zero the 
allowed region is close to the ai-axis. For fx = 100 GeV, about a third of the off-axis 

3 Note that since both fields acquire expectation values, eventually also the second field must undergo 
symmetry breaking, but then as a result of the first field going through its spinodal transition. 




(3.13) 




(3.14) 



poini?. If both eigenvalues are positive, no spinodal instability occurs and we reject the 
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Figure 1: The allowed values of a± 2 after all constraints have been applied, for different 
values of /i. Black dots are maxima, red dots are saddle points. Two superposed points 
refer to the two different allowed values of 0:3. We perform simulations at every second 
allowed point at fx = 100 GeV. 

points are maxima, the rest are saddle points. As fi is further increased, the allowed 
region spreads out to a band near a\ = 0.5 which reconnects around the circle in the 
a2-direction. A "hole" also opens up around the origin. At the largest //, all off-axis points 
are saddle points, and by fi = 400 GeV, no points survive. Interestingly, by far the most 
important constraint is that the minimum should be the global minimum. All but a few of 
the discarded points fail in this respect. We expect that a similar picture arises for other 
choices of 771^2,3, ± with the allowed region shifted accordingly in \i — a\ — a2-space. 

We now turn to our numerical lattice simulations, where we have computed the baryon 
asymmetry for all the allowed parameter space for fi = 100 GeV, top middle of Fig. 1, but 
with a coarser spacing of 7r/20. 

4 Numerical results: 

The action (2.1) is discretized on a lattice and the classical equations of motion derived 
and solved numerically. Starting from a zero-temperature initial condition, we study the 
evolution of the system through the spinodal transition. Observables are averaged over a 
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Figure 2: Left: The Chern-Simons number and Higgs winding numbers in a single initial 
field realization. The inset shows the Higgs field expectation values squared. Right: The 
Chern-Simons number and Higgs winding numbers averaged over the ensemble. 



statistical ensemble of initial realizations, which is by hand C-, P- and CP-symmetric (for 
details, see [1]). The baryon asymmetry is inferred from the anomaly equation (2.4). 

In fact, because we are initially very far from equilibrium, Chern-Simons number is not 
a very clean observable, since it is in general non-integer and can exhibit large oscillations 
at intermediate times. Instead, we consider the Higgs winding number, which coincides 
with Chern-Simons number at late times, is integer throughout and settles much earlier 
into its late-time value. Since we have two Higgs fields, we also have two winding numbers, 
both of which will eventually match Chern-Simons number, and we identify the late-time 
value by the time at which the two agree, irrespective of the value of the Chern-Simons 
number. 

In [1], we studied the dependence of the baryon asymmetry on the strength of C-/P- 
violation. In the present paper, we fix 5c /p = —21 which is close enough to the linear 
regime that we can interpolate to smaller values [1] (^c/P = gives zero asymmetry by 
construction). In this way, we can investigate the significance of CP violation by studying 
the dependence of the baryon asymmetry on 0^2,3 and 8. 

Fig. 2 (left) shows the evolution of winding numbers and Chern-Simons number for one 
particular configuration for a particular choice of (01,02,03) = (0, vr/10, 1.138). Winding 
number has clearly settled, while Chern-Simons number is still catching up. In the inset, we 
show the Higgs field expectation values squared, of which one settles very rapidly, and one 
keeps oscillating for a long time, and with large amplitude. This is because the potential 
around the minimum is steep in one direction and shallow in the other. 

In Fig. 2 (right) we show the average winding number and Chern-Simons number, 
averaged over 100 sets of 4 conjugate configurations. Most of the high-frequency noise in 
the winding numbers has been averaged out, and the two nicely settle at a common value, 
quite early on in the evolution. By vt = 30, symmetry breaking is complete. We also see 
that statistical errors are well under control at this size of ensemble. 

The average Chern-Simons number, however, has certainly not settled to its equilib- 
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Figure 3: Higgs winding number for a parameter set {A} and for {A}*, with an overall 
flipped sign. Here, (ai, «2) 0) = (0, yg, 100 GeV, —1.27) and the results are averaged over 
an ensemble of 4 x 25 configurations. 

rium value. Two effects are at work here: There is a net shift downwards, which is a 
transient non-equilibrium effect. We checked, by running for three times as long, that 
eventually the Chern-Simons number settles to the winding number value. 

The second effect is a large-amplitude oscillation, and is the result of the C-/P-violating 
driving force being given by the oscillating Higgs field vevs. We see from the figure that 
the oscillation has the same frequency as the Higgs field oscillations, but are shifted by a 
phase. This follows from considering the C-/P-violating term as 

5 C/P <x 5 C/P Im^s] d t N cs , (4.1) 

which holds approximately for almost homogeneous Higgs fields. By partial integration, 
this term can be considered a time dependent driving force or chemical potential for Chern- 
Simons number, with magnitude oc <5c/p3t(Im[</>j</>2]). The reason why this second effect is 
not washed out by the averaging procedure is that all members of the ensemble experience 
(roughly) the same oscillation frequency and phase of the driving force, since the Higgs 
oscillation is almost universal, configuration by configuration. Therefore, although other 
configuration-specific effects average out to give a small baryon asymmetry, the driven 
oscillation is common to all configuration and survives the averaging process. At late 
times, the Higgs fields will also stop oscillating, and the driving force will disappear. But 
even at these early times, the coherent oscillation has no impact on the average winding 
numbers, which we therefore take as our measurement of the generated baryon asymmetry. 

4.1 Symmetry under C and A — > A* 

In Fig. 3, we demonstrate explicitly that the symmetry (3.10) holds, by simply computing 
the asymmetry for a parameter set {A} and its complex conjugate, and then flipping the 
sign of the resulting asymmetry. We see that the agreement is very good (within statistical 
error bars). Hence we find the advertised redundancy between positive and negative values 
of Ot2- 
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Figure 4: The dependence of the winding number asymmetry on 9 for (a\,oi2,^) = 
(0, j^, 100 GeV) (black and red) and {?f, fj, 100 GeV) (green and blue). 

4.2 Dependence on 9 

The transformation (3.11) allows us, from a given set of parameters {A}, to generate a 
whole set of identical potentials, but where the minimum is rotated to v\e~ % ® for any value 
of 9. We generate {A} at 9 = using the constraint Im /x^ 2 = wi^ImAs, and from the point 
of view of CP-violation, all values of 9 are equivalent. But once we couple to C-/P-violation, 
the potentials are distinct. In the vacuum, we have 

£ C/P oc 5 C/P 2 vi v 2 sm(9) Tr F^F^, (4.2) 

and so were we in vacuum throughout the transition, there would be no asymmetry gen- 
erated at 9 = 0. And naively, one would expect the asymmetry to be proportional to 
sin#. 

Fig. 4 shows the asymmetry in JY^ as a function of 9 for the parameter sets {A}(ai, a?,, //) 
(0, ^, 100 GeV) (black dots) and {?f, |j, 100 GeV) (green dots). The arrows indicate the 
values of 9 corresponding to real A5, the criterion we will use for most of our simulations 
below. We have fit with a form Asin(9 + 59), and find beautiful agreement, but with a 
non-zero 59 = 2.26 (red line) and 59 = 1.8 (blue line). As a result, even at 9 = 0, an 
asymmetry is generated during the transition where 9 is different from its vacuum value. 

This is a result of the Higgs fields rolling down the potential in a spinodal transition, 
where both the length and the phase of the fields vary locally, until they finally settle 
near their vacuum values. The surprising result is perhaps that the simple sin# form is 
preserved, and that the out-of-equilibrium stage is encoded in the ^-dependence as an 
overall shift of the phase, 59. 

But this also means that the asymmetry vanishes at 9 = —59 (and 9 = ir — 59), and 
that the overall sign of the asymmetry varies in this simple way with 9, presumable for 
any set {A}. We do not know of an obvious way of determining 59 apart from through the 
simulations. On the other hand, since we can parameterize the dependence through A and 
59, we only need simulations at two values of 9, say 9 = and 9 = tt/2. Then we simply 
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Figure 5: The final baryon asymmetry as a function of a\ for ai = ir/20 (left) and n/10 
(right). 

have 

ton(J9) = = ") , A = ^-») , (4.3) 

v ; N w (6 = ir/2) sm(56) y ' 

from which one can find 56 and then A 

In the following, we have for each {A}[0] rotated to the value of 9 that gives a real A5. 
This choice is arbitrary, and can as we have seen with a comparable amount of additional 
some numerical effort be extended to a complete ^-dependence. 

4.3 Dependence on uxp 

Fig. 5 is the baryon asymmetry as a function of a\ at \i = 100 GeV for ai = 7r/20 (left) 
and 7r/10 (right). We have used the conversion from winding number to baryon number 

3S = 1.2 x 10- (M4M)-' ;< x {N y } . (4 . 4) 

This assumes that the total potential energy is distributed onto all the Standard Model 
degrees of freedom with masses less than m w , giving the reheating temperature and the 
photon number density n 7 . 

As discussed in the previous section, the 02 > and 02 < regions are related by 
Eq. (3.9) Since there are at most two allowed values of 03 at each grid point in the («i, 02) 
plane, we finally need to perform our numerical simulation using 39 sets of parameters. 

At 0J2 = it/ 20 (left) we see a clear peak close to a\ = 0, which gradually decreases 
toward the edges of the allowed parameter region. For large a\ the results are roughly 
consistent with zero. The two values of 03 happen to give very similar results within 
errors. Note that connecting the largest /smallest Q3-results by curves is an arbitrary 
choice to guide the eye. At 0.2 = 7r/10 (right) we first of all observe that the two values of 
do not agree as well, although they are still within a factor of two. The asymmetry is 
larger than for 02 = vr/20, and shows no sign of smoothly going to zero at the edge of the 
parameter range. This may be a result of the coarse resolution in a\. Again, the connecting 
curves are just to guide the eye. We see that the maximum value is again attained near 
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ct\ = 0, although now a peak structure is less clear. We should also keep in mind the 
^-dependence, and that the results shown here could be at any point in the period of the 
sin# behaviour. 

The small remaining parameter range at large 02 (see Fig. 1) gives baryon asymmetries 
of roughly the same size, and are by no means suppressed compared to small 0.2- We should 
also mention that we checked that the magnitude of the asymmetry is not in a simple way 
correlated with the determinant or eigenvalues of the mass matrix M. 2 at the origin, and 
in particular whether we start at a maximum or a saddle point. There is also no simple 
correlation with the phase of n\ 2 or of v\. We did, however find a weak correlation in 
the combined Im(ui)-(V r (0, 0) — V(v±, U2))-plane, suggesting that a deep potential drop and 
large CP-violation gives a large baryon asymmetry. This is perhaps not unexpected, but 
is surprisingly difficult to confirm. Clearly, the complicated non-linear dynamics does not 
allow for such simple conclusions about the generated asymmetry. 

5 Conclusion and outlook 

We have outlined a practical parametrization of the 4-dimensional parameter space in the 
2HDM resulting from fixing masses and vevs. We have seen that imposing a number of 
general consistency criteria, a finite region in \x — a\ — «2 survives, and this can be extended 
to the ^-direction by a simple phase change transformation. 

The parameter space is a hyper-cylinder with \i > and three angles a\, 02 and 0, 
with the additional redundancy that 02 > and «2 < are connected. On the other hand 
each set of ot\ t 2 has up to two solutions for 03. When the bosonic sector is coupled to 
fermions, different values of are physically distinct, but the asymmetry seems to follow 
a form Asin(0 + 50). 50 is a priori unknown, but can be found numerically by using two 
different values of 0, say spaced by it/ 2. Finally, it seems that the extent in \i is finite and 
determined by the overall scale of the fixed Higgs masses. 

The maximal asymmetry we found is at the point (01,02,03) = (tt/20, 7t/10, 1.333), 
and is 

(m) = -1.1 x 1(T 5 x 5 C/P , (5.1) 

V n 1 J max 

so that in order to reproduce that observed asymmetry of ~ 6 x 10 -10 , we require 

5 C/P ~ (5 to 6) x 1(T 5 . (5.2) 

We note that the maximal asymmetry is a factor 3 or 4 larger than what we found in [1] at 
one particular parameter point, which therefore was not a particularly unique case. Also, 
there are parts of the allowed parameter space that give vanishing asymmetry. Computing 
5q/p from first principles could therefore potentially rule out regions of 2HDM parameter 
space, under the assumption that baryogenesis originates from a cold spinodal transition 
involving two Higgs field. 

Although a similar programme could be attempted for other scenarios of baryogenesis, 
these are less amenable to a direct, quantitative computation. Leptogenesis is a multi-stage 
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process, generation of lepton asymmetry, thermalization, sphaleron processes, freeze-out. 
And "Hot" electroweak baryogenesis involves the nucleation of bubbles, their interaction 
with the plasma and again sphaleron processes. Cold electroweak baryogenesis offers a 
practicable testing groud for this kind of parameter scans. 

Given the numerical effort involved in the present work (of order 10 5 CPU hours on 
a standard linux cluster), it is difficult to scan through the currently allowed parameter 
space, including the remaining Higgs masses and tan /3. But a complete sweep of the 4- 
dimensional parameter space can be done with about a factor of 10-100 more computing 
power, which is easily within reach of current supercomputers. And hopefully, the coming 
years of LHC-experiments at the electroweak energy scale will constrain the viable range 
of masses and vevs, or even discover additional scalar particles. When this happens, it 
would be natural to revisit the scenario considered here, and use the baryon asymmetry to 
narrow down the range of experimentally less accessible parameters. 
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A Masses of the physical Higgs bosons 

The physical Higgs bosons are defined by the mass eigenstates, which can be found by 
writing 

* = e " ( (-1 + m + ixi) A/2 ) ' 02 = ( (t* + m + i X2 ) /V2 ) " {AA) 

There are in total eight (real) fields but three of them contribute to the degrees of freedom 
of the massive gauge bosons VF ± and Z° after spontaneous symmetry breaking. Let us 
write 

Xi = cos/3G°-sin/3?7 3 , X 2 = cos /V + sin f3G° , (A.2) 
<ft = cos/3G + - sin /3H + , $ = cos/3#+ + sin/3G + , (A.3) 

where G° and (G~ is the complex conjugate of G + ) are Goldstone bosons to be swal- 
lowed up by the gauge bosons. Inserting (A.2) and (A.3) into (2.1), we have for the 
remaining degrees of freedom 

V = —r]M 2 ri T + m±H~ H + + interaction terms, (A. 4) 

where (using sp = sin/3, cp = cos/3, tp = tan/3) 

ml = \ \-v\X, + \f°) + p$*/{c P 8p)] , (A.5) 
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with 



A 5 = X 5 e~ 2W = Af + Uf , and #f 2 = = At<T + *MlT- (A.6) 

are the charged Higgs bosons with mass m-|-. The three neutral physical Higgs bosons 
are obtained by diagonalizing M 2 , which reads explicitly 



/ 2\ ( \ 2 1 1 ~2,Rc, 1 -2,Re . 2\ 



M 



1 » 2, R° i 2 \ 1 * z.ne / , . 2 \ / \ 

/i 12 + t> Z A 345 C / 3S /3 2/^12 IH + V M S /3J 



1 -2,Rc 



V 



1 -2,1m / 
"2^12 l C fi 



1 -2,1m / 
■3^12 



1 -2,1m / 
"2^12 / C /3 

1 -,2,1m / 
"2^12 l S P 



(A.7) 



with A345 = (A 3 



A 4 + Xf c ). Following Ref. [18, 19], we introduce the rotational matrix 

* ClC 2 c 2 si s 2 \ 

-C3S1 - cis 2 s 3 cic 3 - S1S2S3 c 2 s 3 , (A.8 
\ -C1C3S2 + s 1 s 3 -C3S1S2 - C1S3 c 2 c 3 J 



Rn 



such that, 



M = R n diag{m 1 ,m 2 ,m 3 } R n , 



(A.9) 



where Sj = sin a, and Cj = cos on with i = 1,2, 3, and the rotational angles ^ > a±, a 2 > — ^ 
and 7r > 03 > 0. There are 6 independent equations in (A.9), and one of them gives a 
constraint between the three rotational angles, which reads 



(m 3 /m 2 y 



sin(2a3) tan(«i + 0) — 2 sin(a2) ((mi/m 2 ) 2 — sin(«3) 2 ) 



sin(2a3) tan(ai + 0) — 2sin(a2)(l — sin(a3) 2 ) 
and to which there are 0, 1 or 2 solutions for a 3 for each pair (a\,a 2 ). 



(A.10) 
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